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Abstract 

Two-dimensional statistically stationary isotropic turbulence with an imposed uniform scalar 
gradient is investigated. Dimensional arguments are presented to predict the inertial range scaling 
of the turbulent scalar flux spectrum in both the inverse cascade range and the enstrophy cascade 
range for small and unity Schmidt numbers. The scaling predictions are checked by direct numerical 
simulations and good agreement is observed. 
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I. INTRODUCTION 



In the present work we consider the spectral distribution of the passive scalar flux in 
two-dimensional incompressible Navier-Stokes turbulence. The scalar flux appears as the 
unclosed quantity in the Reynolds averaged equation for the mean scalar field: separating 
the velocity and passive scalar field into mean and fluctuations, u = u + u' and 6 = 6 + 6', 
the equation for the mean scalar field reads 



86 _ 06__ d*l_ du' 3 6> 
dt 3 dxj 8x 2 a dxj 

where k is the diffusivity of the scalar and • denotes an ensemble average. The last term 



of this equation contains the correlation u'jO', which is called the scalar flux. It is the term 
which represents the influence of the turbulent fluctuations on the mean scalar profile. Since 
it is the unclosed term in the Reynolds averaged equations, it needs to be modeled, e.g. by 
means of an eddy diffusivity. To propose correct models for the scalar flux, understanding 
of the physics of the turbulent flux is needed. For an overview of models for the scalar 

n 

flux, we refer to the book by Schiestel pj, the work by Rogers 



model derived by Wikstrom et al. [3]. For the more complicated case of the scalar flux in 
the presence of shear and rotation, see the work by Brethouwer [4j]. These studies focus on 
three-dimensional turbulence. 

We consider statistically homogeneous velocity and scalar fields so that we can investigate 
the scale-distribution of the turbulent scalar flux by means of Fourier spectra. The Fourier 
spectrum related to the scalar flux is defined as 

F Uj9 (k)= [ JV^(x,t)0'(^t)]dE(fc), (2) 

in which is a circular wavenumber shell with radius k, the wavenumber, and .F|a;-a:'[-] 
denotes the Fourier transform with respect to the separation vector x — x' . This definition 
is such that by construction we have 

roa 

/ F Ujd (k)dk = ^6 J , (3) 
Jo 

which illustrates that the scalar flux spectrum characterizes the contribution of different 
lengthscales (or wavenumbers) to the scalar flux. This spectrum is also called the scalar- 
velocity co-spectrum since it is defined as the real part of the scalar-velocity correlation 



2| or more recently the 



in Fourier space. The imaginary part is called the quadrature spectrum. The quadrature 
spectrum does not contribute to the scalar flux in physical space and we therefore concentrate 
on the co-spectrum. 

Academically the least complicated case to study the turbulent scalar flux is, as proposed 
by Corrsin [5], isotropic turbulence on which we impose a stationary uniform mean scalar 
gradient d0/dx\ = T, arbitrarily chosen in the Xi-direction. In this case there exists one 
non-zero component of the scalar flux, aligned with the gradient. The other component 
is zero. We consider this case and in particular, we focus on the inertial range scaling of 
the scalar flux spectrum. We will in the following drop the subscripts and denote the co- 
spectrum by F(k). We will also drop the primes and denote the fluctuations of velocity and 
scalar by u and 0, respectively. Before starting the study of the scaling in two-dimensional 
turbulence, we briefly discuss the results obtained in the related case of three-dimensional 
turM eace. L umley 6Q predlcted that at high Re ynolds n Umb e re tbe i^ial r aa g e should 
fall off as k~ 7 / 3 . Indeed he predicted the inertial range to be given by 

F(k) ~ re 1/3 k~ 7 / 3 , (4) 

with e the dissipation of kinetic energy, or more precisely the energy flux at scale k. This 
scaling was investigated experimentally in the atmospheric boundary layer [a| and in de- 
caying grid turbulence at Taylor-scale Reynolds numbers up to R\ = 600 jlllloT]. In these 

re —7/3 scaling was not observed at this 



11] that the inertial range exponent might 



grid-turbulence experiments it was found that t 
Reynolds number. It was subsequently proposed 
be —2 instead of —7/3. However, in closure calculations it was shown that the —2 scaling 
was a low- Reynolds number effect and that the —7/3 scaling should be observed at higher 



Reynolds numbers 
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131 ] . This was confirmed by the work of O'Gorman and Pullin 



and recent direct numerical simulations (DNS) 

In the case of two-dimensional turbulence only few studies address the problem of the 
scaling of the scalar flux spectrum. Let us recall that in two-dimensional turbulence, in 
which the energy is injected at a wavenumber ki, two cascades can be observed, first an 
energy cascade towards the large scales and, secondly, an enstrophy cascade to the small 
scales. If the injection scale is much smaller than the domain-size and much larger than the 
range in which the viscous stresses become important, both cascades are characterized by 
power-law scaling 16Nl8| . We focus on these inertial ranges, which we will denote by IC for 
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the inverse energy cascade and FC for the forward enstrophy cascade range. In particular 
we investigate the wavenumber dependence of the scalar flux spectrum in these ranges. 
One of the few works investigating the scaling of the scalar flux spectrum in two- 



19j . which mentions that the scalar flux spectrum can be roughly 



dimensional turbulence is 
estimated by 

F(k) » Eik^Eeik) 1 / 2 , (5) 
in which the scalar variance spectrum is defined as 

E e {k,t) = \f T\ x ^ x/ [6(x,t)e(x',t)}dJ:(k). (6) 

1 JT,{k) 

In the inverse cascade range, where both the energy spectrum E(k) and the scalar variance 



3,Q 



spectrum E e (k) are known to obey Kolmogorov-Obukhov scaling [20|, |21j], this would lead 
to a fc -5 / 3 inertial range. Close observation of the numerical results in reference [l9| shows 
that this is not the case. 

In the present paper we show that this k~ 5 ^ 3 inertial range prediction does not correspond 
to the physics of the problem. Phenomenological scalings for the inertial ranges in both the 
inverse cascade and the forward enstrophy cascade will be proposed for the scalar flux 
spectrum F(k) and the scalar variance spectrum Eg(k) for the cases of unity and small 
Schmidt number (the Schmidt number is defined as the ratio of the diffusivity of momentum 
to that of the scalar, Sc = v/k, and is identical to the Prandtl number when the passive 
scalar is temperature). Direct numerical simulations are carried out to verify the validity of 
the predictions. 

Note that since the scalar fluctuations are produced by a mean gradient, the scalar 
fluctuations are in principle not isotropic, but axisymmetric around the direction of the 



gradient. It was shown 22|, |23| that in the case of three-dimensional isotropic turbulence 



the spectral distribution of scalar flux can be described by a single scalar function. The 
distribution of scalar variance can be described by two scalar functions. In the present 
work, by integrating over wavenumber shells [equations (E]), (JS}], we eliminate the angle 
dependence. A detailed study of the anisotropy of the scalar field will not be performed in 
the present work. 
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II. LINK BETWEEN THE LAGRANGIAN TIMESCALE AND SCALAR FLUX 
SPECTRUM 



The phenomenological scaling for the scalar flux proposed in the present work is based 
on the direct relation which exists between the scalar field and the Lagrangian dynamics 
of the turbulent velocity field. We therefore first discuss this link. Kraichnan proposed in 
the framework of the Lagrangian History Direct Interaction Approximation 24|, that the 
dominant spectral timescale characterizing the inertial range dynamics can be estimated by 

r(M)= /* ~~Tn7r^~V~ ds — —pr — r f E(k,t\s)ds. (7) 



E(k,t) E(k,t) 
This quantity was investigated numerically in 



o 



251 ] . The energy spectrum is the spherically 



averaged Fourier transform of the two-point velocity correlation, 

E(k,t) = l [ F\ a -*[u i (x,t)ui(x',t)]dE{k). (8) 

E(k,t\s) is the equivalent spectrum in which the Eulerian velocity Ui(x',t) is replaced by 
Ui(x', t\s), which is defined as the velocity at time s of a fluid particle which arrives at point 
x' at time t. The definition of Ui(x', t\s) is illustrated in Figured] The definition of E(k, t\s) 
is thus 

E(k,t\s) = - J F\ x - x i[ui(x,t)ui(x',t\s)]dE(k). (9) 

1 JT.{k) 

By definition E(k, t\t) coincides with the Eulerian spectrum E(k, t). An interesting property 
of (J71) is that the integral can be explicited by integrating Ui(x,t\s) along its trajectory. 

ft 

E(k,t\s)ds = 



\ \ F\x-x'[ui(x,t) I Ui(x',t\s)ds]dZ(k) 

* JS(fc) Jo 

= \ ( J r U_ (B '[ui(aj,*)X i (aj / ,*)]dE(fe). (10) 

1 JUk) 



'E(fc) 

Instead of the two-time quantity Ui(x',t\s), the expression now contains the single-time 
displacement vector of the fluid particle, Xi(x',t), corresponding to the vector pointing 
from its position at t — to its position at t, x', or, in other words, the trajectory. The 
link between the scalar flux spectrum and the integral of E(k,t\s) becomes evident if we 
compare the evolution equation of a non-diffusive passive scalar fluctuation 9 in the presence 
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of a mean scalar gradient, 

dO dO 

with the equation of the xi-component of the Lagrangian position vector Xi(x,t): 

dx 1 ax 1 dx 1 

^r = ^r +u ^ = Ul - (12) 

Indeed, both equations are identical, only differing by a factor —T. As was already stated 



in 26[, the scalar fluctuation is therefore proportional to the displacement of a fluid particle 
in the direction of the gradient. In the limit of vanishing diffusivity, relation ((7J) can thus 
be recasted, using (j2J) and ( TTUl) as in 



3: 



If the energy spectrum and the Lagrangian timescale are known, the scalar flux spectrum is 
given by relation (fl3l) . 



A. Prediction of the scaling of the scalar flux spectrum at large and unity Schmidt 
number 



Dimensional analysis and phenomenological reasoning 24|, |28) give that at a scale I ~ k~ x 
the Lagrangian timescale should be approximately given by l/u(l) in which the typical 
velocity u{l) can be estimated to be of order yjk E{k). This yields an estimation for the 
timescale r(k), 

r(fc) ~ (k 3 E(k)y 1/2 . (14) 

Combining this relation with ( fT3l) yields an estimation for the scalar flux inertial range 
scaling, 

m ~ r^M, (is) 

which is a direct relation between inertial range scaling of the scalar flux spectrum and the 
energy spectrum. In three-dimensional turbulence, using Kolmogorov scaling for the energy 
spectrum, 

E{k) ~ e 2/3 k' 5 / 3 , (16) 
leads to classical scaling for the scalar flux-spectrum, 

F(k) ~ re 1 / 3 *;- 7 / 3 . (17) 



In two-dimensional turbulence this scaling should hold in the inverse cascade range where 
Kolmogorov scaling is expected. In the forward enstrophy cascade range, the energy spec- 
trum is predicted to scale as [16l-ll8|. 



E(k) ~ /3 2/3 A- 3 , 



(18) 



with the flux of enstrophy in the direct cascade. This scaling was later refined introducing 



logarithmic corrections 29 



E(k) ~^ 2/3 fc- 3 /ln(jfc/jfc i ) 1 / 3 , 



(19) 



with hi the wavenumber corresponding to the energy injection. We neglect this correction 
as a first approach. For this forward entrophy cascade range (IT5|) yields the scaling 



F(k) ~ T(3 1/3 k~ 3 . 



(20) 



It should be noted that the preceding analysis supposes a high Schmidt number. Indeed, the 
analogy between the position of a fluid particle and a scalar fluctuation [equations ( ITT]) and 
012p ] is exact for infinite Schmidt number. However, the effect of the Schmidt number for 



Sc larger than one is small 



22 



14| showed that when changing 



31| . O' Gorman and Pullin 
the Schmidt number from 1 to 10 4 , the shape of the scalar flux spectrum was only little 
affected. We now explain this. 

The equation for the co-spectrum can be derived directly from the scalar advection- 
diffusion equation combined with the Navier-Stokes equation (e.g. 
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22j). It reads 



d_ 
dt 



+ (v + K)k? 



F(k) = --YE(k)+T^ L (k). 



(21) 



The left hand side contains the time-derivative and the influence of viscosity v and scalar 
diffusivity k. We consider the statistically stationary state in which the time-derivative term 
drops. The first term on the right hand side is the production of scalar flux by interaction of 
the velocity field with the mean scalar gradient T. The last term is the nonlinear interaction 
which contains two contributions: a purely conservative nonlinear interaction which sums 
to zero by integration over wavenumbers and a purely destructive pressure scrambling term 
which annihilates the correlation between scalar and velocity fluctuations. The viscous- 
diffusive term can be written as 



(u + n)k 2 F(k) = u(l + Sc" )k F(k) 



(22) 



This term changes only by a factor 2 when the Schmidt number goes from 1 to oo. The 
influence of the Schmidt number for Sc larger than one is therefore small. 



B. Prediction of the scaling of the scalar flux spectrum at small Schmidt number 



In the case of Sc — > we do expect the above reasoning to change. We now discuss this 
case of small Schmidt number. 

When the diffusivity becomes very large (keeping v constant to retain an inertial range 
for the energy spectrum), the influence of the nonlinear terms in equation (}2~Tj) will become 
small, since the diffusive timescale becomes smaller than the nonlinear timescale (such as 
the eddy turnover time). The production term is then directly balanced by the diffusive 
term. In this case ( 12T]) reduces to the equilibrium 



which yields 



nk 2 F(k) = —-TE(k) 
2TE(k) 



F{k) 



3nk 2 



(23) 



(24) 



0' Gorman and Pullin [l^] obtained the same expression in three-dimensions. In the inverse 
cascade (IC) this should yield a k~ n ^ 3 scaling, in the forward enstrophy range (FC), a k~ 5 
scaling. 

In section IIVI results of direct numerical simulations of isotropic 2D turbulence with an 
imposed mean scalar gradient are presented to check the relations: 



T- 1 F(A;) ~ < 



e l/3 k -7/3 IC for Sc>l, 

0i/3 k -3 FC for Sc>l, 

^-1^/3^-11/3 JQ for Sc ^ 1; 

K- l (5 2 ^k~ b FC for Sc < 1. 



(25) 



III. PREDICTIONS FOR THE SPECTRUM OF THE PASSIVE SCALAR VARI- 
ANCE 



It is expected that the scalar variance spectrum displays Batchelor scaling 32] in the 
forward enstrophy cascade as was experimentally demonstrated by 33], 



E e (k) ~ e e /r 1/3 £r\ 



(26) 



with eg the (diffusive) destruction rate of passive scalar fluctuations. In the inverse cascade, 
Corrsin-Obukhov scaling is expected. 



Eg(k) ~ tet-^k-^. 



(27) 



The equation for the scalar variance spectrum reads 

E e (k) = -F{k)T + T NL (k), 



j| + 2Kk2 
at 



(28) 



with T e NL (k) the nonlinear transfer term. For very small Schmidt number this equation can 
again be linearized, yielding for the statistically stationary state 

-F(k)T 



This gives, using (J2 



Eg(k) 



E e {k) 



2nk 2 



E(k)T 2 



3K 2 k* 

For the scalar variance, our predictions are therefore 



(29) 



(30) 



Eg(k) ~ { 



e e e- l ^k-^ 3 IC for Sc>l, 
eop-Wk- 1 FC for Sc>l, 
T 2 K- 2 e 2 / 3 k- 17 / 3 IC for Sc < 1, 
k T 2 K- 2 f3 2 /*k~ 7 FC for Sc < 1. 



(31) 



IV. NUMERICAL VERIFICATION OF THE PROPOSED INERTIAL RANGE 
SCALINGS 



Numerical method 



Simulations are performed using a standard pseudo-spectral method [34[]. The simulations 
are fully dealiased and the resolution is 1024 2 gridpoints for a square periodic domain of size 
2/T. The time is advanced using a second order Adams-Bashforth time-stepping scheme. 

The equations for the vorticity field and scalar field are 



du du 



09 
di 



89 

dx~ 



d 2a u d- 2 u 
Va ~dxf + 1 1 faf 

pp.a' a 



(32) 
(33) 



with the vorticity u = e z ■ (V x it), / a random-phase isotropic forcing localized in a band 
in wavenumber-space with a time-correlation equal to the timestep. The parameters a and 
a' are integers equal to one in the case of Newtonian viscosity and diffusivity and equal to 
8 in the case of hyperviscosity or hyperdiffusivity. The mean gradient T is in all cases taken 
equal to 1 so that the scalar flux, and its spectrum, are dominantly negative. 

In all cases, hyperviscosity is used to concentrate the influence of the viscous term at the 
highest wavenumbers. This allows to increase the extend of the inertial range, which is the 
main subject in the present work. Equivalently the scalar variance is removed at the largest 
wavenumbers by a hyper- diffusive term except in the case of small Schmidt number. Since 
in that case the diffusive term becomes the dominant mechanism, the scaling is directly 
affected by the type of diffusion, as can be seen in expressions fl23J) and f[2Uj) . In that case we 
therefore use a 'normal' Laplacian diffusive term (a' = 1). In two-dimensional turbulence 
the energy shows a tendency to cascade to smaller wavenumbers, i.e., to larger scales. To 
avoid a pile-up of energy at the smallest wavenumber linear Rayleigh friction (the last term 
in eq. (|32l) ) is used, with 7 equal to unity. 

Two-different fully developed turbulent flows are investigated. First the inverse cascade 
range, in which the forcing is localized in a wavenumber shell around k; t = 210. In this case 
the forward enstrophy range is reduced to less than an octave and a full decade of inverse 
cascade inertial range is observed in the simulations. Second the forward enstrophy range. 
In this case the forcing is localized around hi = 8, and the inverse cascade range is absent 
since the friction acts strongly in the region k < k{. Parameters used in the simulations 
are summarized in table HI Also shown are some average values of some typical turbulence 
quantities. 

In both velocity fields two different cases are considered for the passive scalar. One at 
Sc = 1, with hyperdiffusivity (a' = 8) and one at small Schmidt number and a' = 1. It is 
not straightforward to define Schmidt numbers for these cases. The precise definition of the 
Schmidt numbers is however not important for the present study, but what is important, is 
the location of the inertial ranges and the ranges where diffusivity becomes important. These 
ranges can be determined as follows. We define a wavenumber fc* at which the nonlinear 
timescale r(k) becomes of the order of the diffusive timescale (nk 2 ) -1 . If k* is in the inertial 
range, we can estimate its value by using expression (fT4j) and the inertial range scalings (j!6p . 
(US]). This yields K ~ (e/zt 3 ) 1 / 4 in the inverse cascade and k* ~ (/3/k 3 ) 1 / 6 in the forward 
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cascade. The wavenumber marks the crossover between an inertial-convective range and 
an inertial-diffusive range. We will call unity Schmidt number cases, these cases in which 
both viscosity and diffusivity mainly act in the last two octaves of the energy and scalar 
spectra, i.e. k* is of the order of the viscous wavenumber, (e/z/ 3 ) 1 / 4 . The direct influence of 
the viscosity and diffusivity is then small for wavenumbers smaller than approximately 100. 
In the case of small Sc, a normal diffusive term is used since the scaling depends directly on 
the Laplacian. The diffusivity is here taken large enough for it to act at all scales, including 
the large scales, i.e. k* is of the order of, or smaller than k e , the wavenumber at which the 
energy spectrum peaks. 

Simulations are performed until a statistically stationary velocity field is obtained. The 
spectra are subsequently obtained by averaging over a time-interval of approximately 300 
time-units, until a relatively smooth spectrum is obtained. This corresponds to 270 T e for 
the IC-range and 540 T e for the FC-range. The large-scale turnover-time T e is here defined 



B. Results 

In Figure [2] visualizations of various quantities are shown at an arbitrary time. It is 
observed that the vorticity field contains clear vortical structures in the forward cascade. In 
the inverse cascade the vorticity field seems almost structureless. However, closer inspection 
shows small vortical structures. Visualization of the stream-function shows more clearly that 
these structures are present. The scalar field shows how fluctuations of passive scalar are 
created by interaction of the flow with the mean scalar gradient. In the IC case this scalar 
field is almost structureless, but shows patches of scalar fluctuation. We also displayed the 
instantaneous scalar flux, which is the product of the Xi-component of the velocity with the 
scalar field. Both positive and negative values of the flux are observed. The mean value 
is however smaller than zero (since the mean gradient is positive), so that the net flux is 
non-zero. 

In Figure |3] visualizations are shown for the scalar field and the scalar flux for the small 
Schmidt number case. Vorticity fields and stream function are not shown, since they are 
qualitatively the same as in Figure EJ Due to the large diffusivity, all scalar gradients are 
rapidly smoothed out, so that in both the IC and FC case the scalar field consists of large 
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blobs. The scalar flux fields are characterized by a finer structure. 

In Figure H] wavenumber spectra are shown for the energy, scalar variance and scalar 
flux. In the IC case, classical Kolmogorov scaling proportional to k~ 5 ^ 3 holds for E(k) 
in the inertial range. The scalar variance spectrum Eg(k) is also proportional to k~ 5 ^ 3 as 
can be expected from Corrsin-Obukhov arguments, but showing an important pre-diffusive 
bump. This bump is frequently observed in spectra of the scalar variance, e.g. 0, [3^]. The 
scalar flux spectrum is proportional to k~ 7 ^ 3 which is in disagreement with expression (jSJ) 



proposed as a rough estimate by Smith et al. 19], and in perfect agreement with expression 



(|2"5|) . which corresponds to classical Lumley scaling. Zero-crossings are observed so that not 
the whole spectrum has the same sign. 

In the FC range, the energy spectrum is approximately proportional to k~ 3 , but slightly 
steeper for the wavenumbers close to the injection scale fcj. Taking into account the log- 
arithmic correction, the agreement with the prediction improves even more. The scalar 



variance spectrum Eg(k) shows a Batchelor regime 32| proportional to k . The scalar flux 
spectrum does show a scaling close to the scaling of the energy spectrum, especially for the 
absolute value of the spectrum. It is observed that the spectrum changes sign at several 
wavenumbers. These sign-changes were also observed in the investigation of the scalar flux 



by the stretched spiral vortex model for three-dimensional turbulence [23j . The spectrum of 
the planar contribution of the Lundgren vortex to the scalar flux showed equivalent negative 
excursions. We therefore relate this behavior to the roll-up of the scalar field by large co- 
herent vortices. Indeed, a fluid particle which remains for a long time trapped in a vortical 
structure will contribute both positively and negatively to the scalar flux. 

As can be observed in Figure [5j at small Schmidt number, excellent agreement is observed 
with the predictions. In the IC range, F(k) is proportional to fc -11 / 3 and Eg(k) to k~ X7 l 3 . 
In the FC range, F(k) is proportional to k~ 5 and Eg(k) to k~ 7 . 

V. CONCLUSION 

In this work the scaling of the scalar flux spectrum in two-dimensional isotropic turbulence 
was addressed. Phenomenological arguments based on Lagrangian dynamics were proposed 
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leading to the following predictions for the inertial range scaling of the scalar flux spectrum, 



gl/3^-7/3 



T- L F(k) ~ { 



IC for Sc > 1, 

01/3J.-3 p C for Sc>l, 

^-1^/3^-11/3 IC for Sc < X 

^-1^2/3^-5 FCforSc<l, 



(34) 



and for the scalar variance spectrum, 

Eg{k) ~ < 



IC for 5c > 1, 
FC for Sc > 1, 



(35) 



T 2 K~h 2 / 3 k- 17 / 3 IC for 5c < 1, 
k r 2 /t- 2 /3 2 / 3 A;- 7 FCforSc<l. 

It was shown by DNS that in the inverse cascade the scalar flux spectrum is proportional to 
A;~ 7 / 3 , in perfect agreement with the scaling arguments. The scalar variance shows Corrsin- 
Obukhov scaling, proportional to k~ 5 ^ 3 . In the direct enstrophy cascade the energy-spectrum 
obeys a log-corrected k~ 3 scaling and the scalar spectrum displays Batchelor scaling propor- 
tional to k~ l . The scalar flux spectrum shows important positive and negative contributions, 
probably related to the presence of long-living coherent structures. The absolute value of 
the spectrum shows a scaling close to k~ 3 . At small Schmidt number, excellent agreement is 
observed with the predictions. The scalar flux spectrum scales here as k~ ll l 3 in the IC case 
and k~ 5 in the FC case. The scalar spectrum is proportional to k~ 17 ^ 3 (IC) and k~ 7 (FC). 
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FIG. 2: (Color online) Visualizations of (from top to bottom) vorticity, streamfunction, scalar 
fluctuations, scalar flux. Left: inverse cascade. Right: forward cascade. The Schmidt number is 
unity. The mean scalar gradient is in the horizontal direction. 
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FIG. 3: (Color online) Visualizations of scalar fluctuations (top), scalar flux (bottom), in the 
inverse cascade (left) and in the forward cascade (right) for the case of small Schmidt number. 
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FIG. 4: (Color online) The energy spectrum, scalar flux spectrum and scalar variance spectrum for 
Sc = 1. Top: the case of large wavenumber forcing (inverse energy cascade). Bottom: the case of 
small wavenumber forcing (forward entrophy cascade). The solid lines are dimensional predictions 
given by equation (|25p and (|3ip . In the FC case also the log-corrected k~ s scaling is shown for the 
energy spectrum, which almost superposes the normal A: -3 scaling. Dots indicate positive values 
of the scalar flux spectrum. 
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FIG. 5: (Color online) The energy spectrum, scalar flux spectrum and scalar variance spectrum for 
Sc <C 1. Top: the case of large wavenumber forcing (inverse energy cascade). Bottom: the case of 
small wavenumber forcing (direct entrophy cascade). The solid lines are dimensional predictions 
given by equation (|25p and (|3ip . In the FC case all predictions are also shown with logarithmic 
corrections, which almost superpose on the uncorrected scalings. 
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